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Abstract 

It is speculated how dark energy in a brane world can help reconcile an infinitely 
cyclic cosmology with the second law of thermodynamics. A cyclic model is de- 
scribed, in which dark energy with w < — 1 equation of state leads to a turnaround 
at a time, extremely shortly before the would-be Big Rip, at which both volume 
and entropy of our universe decrease by a gigantic factor, while very many indepen- 
dent similarly small contracting universes are spawned. The entropy of our model 
decreases almost to zero at turnaround but increases for the remainder of the cycle 
by a vanishingly small amount during contraction, empty of matter, then by a large 
factor during inflationary expansion. 



One of the oldest questions in theoretical cosmology is whether an infinitely oscillatory 
universe which avoids an initial singularity can be consistently constructed. As realized by 
Friedmann [1] and especially by Tolman [2,3] one principal obstacle is the second law of 
thermodynamics which dictates that the entropy increases from cycle to cycle. If the cycles 
thereby become longer, extrapolation into the past will lead back to an initial singularity 
again, thus removing the motivation to consider an oscillatory universe in the first place. 
This led to the abandonment of the oscillatory universe by the majority of workers. 

Nevertheless, an infinitely oscillatory universe is a very attractive alternative to the 
Big Bang. One new ingredient in the cosmic make-up is the dark energy discovered only 
in 1998 and so it natural to ask whether this can avoid the difficulties with entropy. 

Some work has been started to exploit the dark energy in allowing cyclicity possibly 
without the need for inflation in [4-8]. Another new ingredient is the use of branes and 
a fourth spatial dimension as in [9-12] which examined consequences for cosmology. The 
Big Rip and replacement of dark energy by modified gravity were explored in [13, 14]. 

If the dark energy has a super- negative equation of state, u\ = p\/ p\ < — 1, it leads 
to a Big Rip [15] at a finite future time where there exist extraordinary conditions with 
regard to density and causality as one approaches the Rip. In the present article we explore 
whether these exceptional conditions can assist in providing an infinitely cyclic model. 

We consider a model where, as we approach the Rip, expansion stops due to a brane 
contribution just short of the Big Rip and there is a turnaround at time t — tr when 
the scale factor is deflated to a very tiny fraction (/) of itself and only one causal patch 
is retained, while the other l// 3 patches contract independently to separate universes. 
Turnaround takes place an extremely short time (< 10~ 27 s) before the Big Rip would 
have occurred, at a time when the universe is fractionated into many independent causal 
patches [14]. 

We discuss contraction which occurs with a very much smaller universe than in expan- 
sion and with almost vanishing entropy because it is assumed empty of dust, matter and 
black holes. 

A bounce takes place a short time before a would-be Big Bang. After the bounce, 
entropy is injected by inflation [16], where is assumed that an inflaton field is excited. 
Inflation is thus be a part of the present model which is one distinction from the work of 
[5-8]. For cyclicity of the entropy, S(t) = S(t + r) to be consistent with thermodynamics it 
is necessary that the deflationary decrease by f 3 compensate the entropy increase acquired 
during expansion including the increase during inflation. 

A possible shortcoming of the proposal could have been the persistence of spacetime 
singularities in cyclic cosmologies [17], but to our understanding for the model we outline 
this problem is avoided, provided that the time average of the Hubble parameter during 
expansion is equal in magnitude and opposite in sign to its average during contraction. 

This model is published because it gives renewed hope for the infinitely oscillatory 
universe saught in [1-3]. Time will tell whether the present model is consistent, but at 
present we see no fatal flaw. 
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Friedmann equation for expansion phase. Let the period of the Universe be designated by 
r and the bounce take place at t — and turnaround at t = t^. Thus the expansion phase 
is for times < t < tj- and the contraction phase corresponds to times tx < t < r. We 
employ the following Friedmann equation for the expansion period < t < t^: 



a(t)\ 2 _ 8nG 
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where the scale factor is normalized to a(to) = 1 at the present time t = to — 14Gy. To 
explain the notation, (pi)o denotes the value of the density pi at time t = t . The first two 
terms are the dark energy and total matter (dark plus luminous) satisfying 
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where H = a(to) / a(t ). The third term in the Friedmann equation is the radiation density 
which is now Q r = 1.3 x 10~ 4 . The final term ~ Ptotai{t) 2 is derivable from a brane set- 
up [9,10,12]; we use a negative sign arising from negative brane tension (a negative sign can 
arise also from a second timelike dimension but that gives difficulties with closed timelike 
paths), ptotai = ^i=A,m,rPi- As the turnaround is approached, the only significant terms 
in Eq.fJTJ) are the first (where u\ < —1) and the last. As the bounce is approached, the 
only important terms in Eq.fJTJ) are the third and the last. (We shall later argue that the 
second term, for matter, is absent during contraction.) In particular, the final term of Eq. 
CQ), ~ Ptotai (t) 2 , arising from the brane set-up is insignificant for almost the entire cycle 
but becomes dominant as one approaches t — > tr for the turnaround and again for t — *■ r 
approaching the bounce. 

Turnaround. Let us assume for algebraic simplicity u\ = —4/3 = constant. This value is 
already almost excluded by WM AP3 [18] but to begin we are aiming only at consistency of 
infinite cyclicity. More realistic values may be discussed elsewhere. With the value uj\ = 
—4/3 we learn from [13] that the time to the Big Rip is (t r i P — to) = HGy(— U\ — 1) _1 = 
33Gy which is, as we shall discuss, within 10 -27 second or less, when turnaround occurs at 
t = t T . So if we adopt t = 14Gy then t T = t + (t rip — to) ~ (14+33)Gfy = 47Gy. From the 
analysis in [13-15] the time when a system becomes gravitationally unbound corresponds 
approximately to the time when the dark energy density matches the mean density of the 
bound system. For an object like the Earth or a hydrogen atom water density pH 2 o is a 
practical unit. 

With this in mind, for the simple case of uo = —4/3 we see from Eq.flTJ that the 
dark energy density grows proportional to the scale factor p\{t) oc a(t) and so given that 
the dark energy at present is pa ~ 10~ 29 g/cm 3 it follows that Ph{t H2 o) = Ph 2 o when 
o(tff 2 o) ~ 10 29 . We can estimate the time tu 2 o by taking on the RHS of the Friedmann 
equation only dark energy (-) = H^Vt^a' 13 with (3 = 3(1 + uj). When we specialize to 
lo = —4/3 it follows that 
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so that (t r ip — tn 2 o) — 33Gy x 10 



-14.5 



10 3,5 s ~ 1 hour. [The value is sensitive to u] It 



is instructive to consider approach to the Rip a more general critical density p c = i]p H2 o 
and to compute the time (t rip — t v ) such that p\{t v ) = p c = r\pu 2 o- We then find, using 
a(t v ) = 10 29 ?y, that 
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which is the required result. We shall see rj > 10 31 so the time in (jlj) is < 10 s. 

To discuss the turnaround analytically we keep only the first and last terms, the only 
significant ones, on the RHS of Eq.(pQ) which becomes for the special case uj = 4/3 
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Writing a = z 2 and z = (oti/ l a 2 ) 1 ^ 2 sin9 gives 

dt = -S£ r— = 
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Integration then gives for the scale factor 
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At turnaround t = t T , a(t T ) = [p c /(px)o] = {a{t)) m ax- At the 
1 and sin 2 9 = [(/?a)q/pc] 1, increasing during subsequent 



where C = -(3/2ttG Pc )^ 2 
present time t = t , a(t ) - 
expansion to 6t = 7r/4. 

A key ingredient in our model is that at turnaround t = tr our universe deflates dra- 
matically with efffective scale factor a(tx) shrinking before contraction to a(tr) — / q (^t) 
where / < 10 -28 . This jettisoning of almost all, a fraction (1 — /), of the accumulated 
entropy is permitted by the exceptional causal structure of the universe. We shall see 
later that the parameter rj at turnaround lies in the range r\ = 10 31 to n — 10 87 which 
implies a dark energy density at turnaround (Planckian density of pa ~ 10 104 pH 2 o can 
be avoided) such that, according to the Big Rip analysis of [13,14], all known, and yet 
unknown smaller, bound systems have become unbound and the constituents causally dis- 
connected. Recall that the density of a hydrogen atom is approximately pH 2 o an d we are 
reaching a dark energy density of from 31 to 87 orders of magnitude higher. 

According to these estimates, at t = tx the universe has already fragmented into an 
astronomical number (l// 3 ) of causal patches, each of which independently contracts as a 
separate universe leading to an infinite multiverse. The entropy at t = tx is thus divided 
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between these new contracting universes and our universe retains only a fraction f 3 . Since 
our model universe has cycled an infinite number of times, the number of parallel universes 
is infinite. 

Friedmann equation for contraction phase. The contraction phase for our universe occurs 
for the period tx < t < r. The scale factor for the contraction phase will be denoted by 
a(t) while we use always the same linear time t subject to the periodicity t + r = t. At 
the turnaround we retain a fraction f 3 of the entropy with a(tx) = fa(tr) and for the 
contraction phase the Friedmann equation is 



h(t) \ 8nG 



a(t) 

where we have defined 
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but in contrast to Eq.([T|) we have set p m = because our hypothesis is that the causal patch 
retained in the model contains only dark energy and radiation but no matter including no 
black holes. This is necessary because during a contracting phase dust or matter would 
clump, even more readily than during expansion, and inevitably interfere with cyclicity. 
Perhaps more importantly, presence of dust or matter would require that our universe 
go in reverse through several phase transitions (recombination, QCD and electroweak to 
name a few) which would violate the second law of thermodynamics. We thus require that 
our universe comes back empty! Any tiny entropy associated with radiation is constant 
during adiabatic contraction. 

The contraction of our universe will proceed from one of the 1/f 3 causal patches 
following Eq.fjH]) until the radiation balances the brane tension at the bounce. 

Bounce. At the bounce, the contraction scale is given, using p c = t]Ph 2 o, from Eq. (JTJ as 

= ( u ) 

Now the model's bounce at t — r must be before the electroweak transition at tsw — 
lCT 10 s when a(tEw) = 1CT 15 , and after the Planck scale when a(tpi anc k) = 10~ 32 in order 
to accommodate the well established weak transition and to avoid uncertainties associated 
with quantum gravity. With this in mind, here are three illustrative values (A, B, C) for 
the bounce temperature T B : 

• (A) At a GUT scale T B = 10 17 GeV,a(t B ) = 1(T 30 . 

• (B) At an intermediate scale T B = 10 10 GeV, a(t B ) = 10" 23 . 

• (C) At a weak scale T B = 10 3 GeV, a(t B ) = 1(T 16 . 
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From Eq. (fTTj) and Eq.Q for these three cases one finds 
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Immediately after the bounce, we assume that an inflaton field is excited and there is 
conventional inflation with enhancement E = a(r + 5)/a(r). Successful inflation requires 
E > 10 28 . Consistency requires therefore / < E^ 1 to allow for the entropy accrued during 
expansion after inflation. The fraction of entropy jettisoned from our universe at deflation 
is thus extremely close to one, being less than one and more than (1 — 10 -28 ) 3 . 

Entropy. Consider first the present epoch t — to- The contributions of the radiation to 
the entropy density s follows the relation 

2vr 2 

8 = —g,T 3 (12) 
45 y K J 

Photons contribute g* = 2he present CMB temperature is T = 2.73K = 0.235meV ~ 
1.191(mm) -1 . Substitution in Eq. (|T2|) gives a present radiation entropy density s 7 (t ) — 
1.48(mm) -3 . Using a volume estimate V = (4it/3)R 3 with R = OGly ~ xl0 29 mm gives a 
total radiation entropy S 7 ~ 6.3 x 10 87 . Including neutrinos increase g* in Eq. fll2p from 
g* = 2 to = 3.36 = 2 + 6 x (7/8) x (4/11) 4 / 3 . This increases S 7 = 6.3 x 10 87 to 
S y+ , ~ xlO 88 . 

This total entropy is interpretable as exp(10 88 ) degrees of freedom, or in information 
theory [19] to a number / of qubits where 2 1 = e s so that I = S/ (ln2 = 0.693) ~ 10 88 . 
This is well below the holographic bound which is dictated by the area in terms of Planck 
units 10~ 64 mm 2 which gives Shoi og (to) = 47r(10 29 mm) 2 /(10~ 32 mm) 2 ~ 10 123 about 10 35 
times bigger. In [19] it is suggested that at least some of this difference may come from 
supermassive black holes. The entropy contribution from the baryons is smaller than S 1 
by some ten orders of magnitude, so like that of the dark matter, is negligible. 

What is the entropy of the dark energy? If it is perfectly homogeneous and non- 
interacting it has zero temperature and entropy. Finally, the 4th term in Eq.([T]) corre- 
sponding to the brane term is neglible, as we have already estimated. The conclusion is 
that at present S to tai{to) ~ 10 88 . 

Now consider the entropy approoaching turnaround at t — tr- We have estimated 
that a{t T ) = 10 29 r] and representative values for r\ = p c / Ph 2 o are 10 31 , 10 59 and 10 87 . 
The temperature T 7 of the radiation scales as T 7 oc a{t)~ l so using the entropy density 
of Eq. (|T2|) a comoving 3- volume oc a(t) 3 will contain the same total radiation entropy 
S-yitr) = S y (to) as at present; this is simply the usual adiabatic expansion. The expansion 
from t = to tx is not purely adiabatic because irreversible processes take place. The 
first is inflation which increases entropy by > 10 84 . There are phase transitions such as 
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the electroweak transition at tsw ~ lOOps, the QCD phase transition at tQCD ~ 100//S, 
and recombination at t rec ~ 10 13 s. Further irreversible processes occur during during 
stellar evolution. Although the expansion of the radiation, the dominant contributor to 
the entropy, is adiabatic, the entropy of matter increases in accord with the second law 
of thermodynamics. In our model, the entropy of the matter increases between t = 
and t T ~ 47Gy. Setting the entropy of the dark energy to zero and the radiation as 
adiabatic, the matter part represented by p m will cause the entropy to rise from S(t = 0) 
to S(t T ) = S(t = 0) + AS where AS causes the contradiction plaguing previous oscillatory 
model universes [1-3]. 

Our main point is that in order for entropy to be cyclic, the entropy which was enhanced 
by a huge factor E 3 > 10 84 at inflation must be reduced dramatically at some point during 
the cycle so that S(t) = S(t + r) becomes possible. Since it increases during expansion 
and contraction, the only logical possibility is the decrease at turnaround as accomplished 
by our causal patch idea. The second law of thermodynamics continues to obtain for 
other causal patches, each with practically vanishing entropy at turnaround, but these are 
permanently removed from our universe contracting instead into separate universes. 

For contraction tx < t < r we are assuming the universe during contraction is empty of 
matter until the bounce so its entropy is vanishingly small. Immediately after the bounce 
inflation arises from an inflaton field, assumed to be excited. We find the counterpoise 
of inflation at the bounce and deflation at turnaround an appealing aspect of the present 
model. 

Conclusion. The standard cosmology based on a Big Bang augmented by an inflation- 
ary era is impressively consistent with the detailed data from WMAP3 [18] when dark 
energy, most conservatively a cosmological constant, is included. Our objections to this 
standard model are more aesthetic than motivated directly by observations. The first 
objection is the nature of the initial singularity and the initial conditions. A second ob- 
jection, not of concern to all colleagues, is that the predicted fate of the universe is an 
infinitely long expansion. We have outlined here a cyclic cosmology resting on phantom 
dark energy where these objections are ameliorated: the classical density and temperature 
never become infinite and future expansion is truncated. Also, our proposal of deflation 
naturally leads to a multiverse picture, somewhat reminiscent of that predicted in eternal 
inflation, though here the proliferation of universes must be infinite and originates at the 
opposite end of a cyclic cosmology, at its maximum rather than at its minimum size. 

We publish our infinitely cyclic model mainly in the hope that it will stimulate a more 
detailed and compelling formulation. 
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